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Abstract. STAR's measurements of directed flow (vi) for pions, kaons (charged and K^), 
protons and antiprotons for Au + Au collisions at ^snn = 200 GeV are presented. Negative 
vi{y) slope is observed for pions, antiprotons, protons and kaons. The proton vi{y) slope at 
midrapidity is found extremely small. Sizable difference is seen between vi of protons and 
anti-protons in 5-30% central collisions. Anti-flow can explain the negative slope, however, it 
has difficulties explaining the centrality dependence of the difference between the vi{y) slopes 
of protons and anti-protons. The vi excitation function is presented. Comparisons to model 
calculations (RQMD, uRQMD, AMPT and QGSM) are made, and it is found that none of the 
four modeis can successfully describe the data. 



1. Introduction 

In relativistic heavy ion cohisions, anisotropic flow describes the collective motion of particles 
with respect to the reaction plane, and it is conventionally characterized by Fourier coefficients 



where (f) denotes the azimuthal angle of an outgoing particle, ip is the orientation of the reaction 
plane, and n denotes the harmonic. So far three such coefficients are found finite at RHIC, 
namely, directed flow vi, elliptic flow V2 and 4*'* order harmonic flow V4. This paper will focus 
on the directed flow, the flrst Fourier coefficient. Directed flow describes the sideward motion 
of fragments in ultra-relativistic nuclear collisions, and it carries early information from the 
collision [2]. The shape of directed flow at midrapidity is of special interest because it might 
reveal a signature of a possible phase transition from normal nuclear matter to a Quark-Gluon 
Plasma (QGP) OHIIS]. It is argued that directed flow, as a function of rapidity, may exhibit a 
flatness at midrapidity due to a strong, tilted expansion of the source. Such expansion gives rise 
to anti-flow ^ or a 3'"'^ flow [1] component. The anti-flow (3'''^ flow component) is perpendicular 
to the source surface, and is in the direction that is opposite to the direction of the bouncing-off 
motion of nucleons. If the tilted expansion is strong enough, it can even overcome the bouncing- 
off motion and results in a negative f i(y) slope at midrapidity, producing a wiggle-like structure. 
In 1^, it is emphasized that proton vi can serve as a probe for the first order phase transition, 
as calculation shows that proton vi (y) does not show negative slope for a hadronic Equation of 
State (EoS) without a QGP phase transition. Note that although the calculation is done with a 
first order phase transition at SPS energies, the direct cause of the negative slope is the strong, 




(1) 



tilted expansion, which is also relevant at top RHIC energies. A wiggle structure is also seen in 
RQMD [6], and it is attributed to the baryon stopping and positive space-momentum correlation. 
In this picture, no phase transition is needed, and pions and nucleons flow in opposite directions. 
To distinguish between baryon stopping and anti-flow associated with a phase transition, it is 
desirable to measure the vi{y) for identified particles and compare the signs of their slope at 
midrapidity. The centrality dependence of proton vi{y) may have implications for a possible 
first order phase transition. It is expected that in very peripheral collisions, protons flow in 
the same direction as spectators. In mid-central collisions, if there is a phase transition, the 
proton vi{y) slope at midrapidity may change sign and become negative. Eventually the slope 
diminishes in central collisions due to the symmetry of collisions. 

At low energies, it is shown by the E895 collaboration [7j that has a negative vi{y) 
slope around midapidity, while lambda and protons have positive slope [8]. It is explained 
by a repulsive kaon-nucleon potential and attractive lambda-nucleon potential. The NA49 
collaboration [9j has measured the vi for pions and protons, and negative vi{y) slope is observed 
by the standard event plane method. The three particle correlation method fi{3}, which is 
believed to be less sensitive to non-flow effects, gives negative slope too, but with larger statistical 
error. At top RHIC energies, vi has been studied for charged particles by both the STAR and 
the PHOBOS collaborations \10\ [TT| \T7[ 113] . It is found that vi in the forward region follows 
the limiting fragmentation hypothesis well, and vi{r]) depends only on the incident energy, but 
not on the size of the colliding system at a given centrality. The system size independence 
of vi can be explained by Hydrodynamic calculation with an titled initial condition |14) . The 
systematic study of vi for identified particles at RHIC did not begin until recently because it 
is more challenging for two reasons, 1) vi for some identified particles (for example, protons) is 
much smaller than that of charged particles, thus is more difficult to measure , 2) it demands 
more statistics to measure vi for identified particles other than poins. 

2. Data Set 

In this analysis we used 62 million events for Au -|- Au collisions at y^SNN = 200 GeV, all taken 
by star's minimum bias trigger during RHIC run 2007. Charged particles are identified by 
their energy loss inside STAR's Time Projection Chamber (TPC) [TS], and -fr^(— t- vr'^vr") are 
reconstructed by their charged daughter tracks inside the TPC. Track quality cuts are the same 
as used in [T6]. The centrality definition of an event was based on the number of charged global 
tracks in the TPC with track quality cuts: |ry| <0.5, a Distance of Closest Approach (DCA) to 
the global vertex less than 3 cm, and 10 or more fit points. Additional weight is assigned to 
each event in the analysis accounting for the bias in the vertex Z direction, due to the inclusion 
of hits of Silicon Vertex Tracker in tracking for year 2007 data production. In addition, the 
transverse momentum pT for protons are required to be larger than 400 MeV/c, and DCA are 
required to be less than 1 cm, in order to avoid the beam-pipe background. The same cut is 
applied to anti-protons as well to enssure a fair comparison with protons. The high end of the 
Pt cut is 1 GeV/c where protons and pions have the same energy loss in the TPC and thus 
become indistinguishable. The event plane angle is determined from the sideward deflection of 
spectator neutrons measured by STAR's Shower Maximum Detector at Zero Degree Calorimeters 
(ZDCSMD), which are located close to beam rapidity and have the minimum contribution from 
non-flow. The description of measuring vi using the ZDCSMD event plane can be found in [12j . 

3. Result 

Fig. [T] shows the pion vi as a function of rapidity. The result is compared to four model 
calculations, namely, RQMD [6j, UrQMD [IT], AMPT [IS] and QGSM [19J. Errors in this plot 
(and other data figures in this paper) are statistical only. Systematical uncertainties are studied 
by varying track quality cuts, and are found to be negligible. Calculations from RQMD and 
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Figure 1. Pion vi as a function of rapidity (or rj, for the QGSM model only) for 10-70% Au 
+ Au collisions at ^snn = 200 GeV. Also plotted are corresponding model calculations. The 
calculation from UrQMD and QGSM are scaled by a factor of 2.5 and 0.18, respectively. Errors 
are statistical only. 



QGSM are scaled by a factor of 2.5 and 0.18, respectively. Following convention, the sign of 
spectator vi in the forward region is chosen to be positive, to which the measured sign of vi 
for particles of interest is only relative. The pion vi{y) result follows closely to the published 
charged particle result |12j and it shows a negative slope at midrapidity. Most models either 
predict the wrong sign of pion vi or the wrong magnitude by a signifiant factor 0, which indicates 
that it is a challenging job to understand the dynamics happening at early stages. In this sense 
the measurement presented in this paper will offer strong constraints on models. 

To address the shadowing effect and the anti-flow, in FiglJl proton and anti-proton vi are 
plotted together with pion vi. Pions are the dominant particles produced and they carry the 
flow of the bulk. It is not a surprise to see that anti-protons, also as produced particles, follow 
the flow direction of pions. They both have negative slope at midrapidity, which is consistent 
with the anti-flow picture. vi{y) of protons, however, exhibits a flatter shape than others. The 
mass difference between protons and pions cannot explain the flatness of proton vi{y) because 
the anti-protons have the same mass as that of protons, yet they flow close to pions. Indeed, 
the observed vi for protons is a convolution of directed flow of produced protons with that 
of transported protons, and the flatness of inclusive proton vi{y) around midrapidity can be 
explained by the negative flow of produced protons being compensated by the positive flow 
of protons transported from spectator rapidity, as a feature expected in the anti-flow picture. 
Fig. [3] shows vi of kaon and as a function of rapidity, for centrality 10-70%. Pion vi is 
also plotted for comparison. Kaon vi is a unique measurement because the kaon/proton cross 
section is smaller than the pion/proton cross section, thus its directed flow is less likely to be 
"induced" by the shadowing effect. The plot shows that both kaon and have negative vi 



^ However, in a recent paper [H] posted after the 26"' Winter Workshop on Nuclear Dynamics, it is shown that 
the directed flow of charged particles can be explained by Hydrodynamic calculation with a tilted initial condition. 
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Figure 2. Proton (solid circles) and anti-proton (solid squares) vi as a function of rapidity 
for 10-70% Au + Au collisions at ^snn = 200 GeV, pion vi{y) (empty stars) is shown as a 
reference. 
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Figure 3. Charged kaon (up triangles) and (down triangles) vi as a function of rapidity 
for 10-70% Au + Au collisions at -y/sNN = 200 GeV, pion vi{y) (empty stars) is shown as a 
reference. 



slope, which means that the observed vi is not likely to be caused by the shadowing effect. 

On the other hand, anti-flow has difficulties in explaining the centrality dependence of vi, 
as shown by Fig. S] , in which vi slope at midrapidity is plotted as a function of centrality for 
protons, anti-protons and charged particles. If there is an anti-flow due to the strong, tilted 
expansion, one expects such an effect is larger in mid-central collisions than that in peripheral 
collisions. As a consequence, proton vi slope, which is expected to be positive in very peripheral 
collisions, will change its sign to negative in mid-central collisions and approach zero in central 
collisions. Due to the large statistical error, we cannot exclude the possibility that proton vl 
can be positive in 70-80% central collisions. However, in 30-80% central collisions, proton vi 
slope is found mostly negative and the magnitude decreases with decreasing centrality. In more 
central (5-30%) collisions, proton vi slope becomes extremely small, while anti-proton vi slope 
remains negative and continues to follow that of charged particles (mostly pions). Anti-flow may 
cause a difference between vi of protons and anti-protons, and such difference is expected to 
be accompanied by strongly negative vi. In data, the large difference between proton and anti- 
proton Vl is seen in 5-30% centrality, while strongly negative vi is found for protons, anti-protons 
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Figure 4. Charged (solid stars), proton (solid circles) and anti-proton (solid squares) vi{y) 
slope (dvi/dy) at midrapidity as a function of centrality for Au + Au collisions at ^snn = 200 
GeV. 



and charged particles in a different centrality (30-70%). 
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Figure 5. Proton vi(y') slope {dvi/dy) at midrapidity as a function of center of mass collision 
energy, where y = y/ybeam- 

In Fig. O proton vi{y / ybeam) slope at midrapidity is plotted as a function of collision energy. 
We see that proton vi slope decreases rapidly with increasing energy, reaching zero around 
y^SNN = 9 GeV. It changes its sign to negative as shown by the data point at -^snn = 17 GeV, 
measured by NA49 experiment. With previous measurements, which include only one point 
above -^/snn = 9 GeV, one cannot conclude if proton vi slope continues to decrease or stays 



close to zero when energy increases, the addition of the data point from RHIC indicates that 
proton f 1 slope remains close to zero at -^/snn = 200 GeV. Judging over the broad energy range, 
the transition of the trend happens around y^SNN = 9 GeV, making it a discontinuous point, 
and interestingly, it coincides with the energy vicinity where {k~^)/{'K^) exhibits a kink [20j. 
Hopefully the on-going Beam Energy Scan program will map out the region of interest with 
details. 

4. Summary 

For 10-70% central Au -|- Au collisions at ^snn = 200 GeV, directed flow for pions, kaons 
(charged and K^), and anti-protons are found to have negative slope at mid-rapidity. Protons, 
however, exhibit a flat shape around mid-rapidity. Sizable difference is seen between vi of 
protons and anti-protons in 5-30% central collisions. Anti-flow can explain the negative slope, 
however, it has difficulties explaining the centrality dependence of the difference between the 
vi{y) slopes of protons and anti-protons. 
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